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Abstract 

The article presents calculated dissociative recombination (DR) rate coefficients for Hg . The 
previous theoretical work on was performed using the adiabatic hyperspherical approximation 
to calculate the target ion vibrational states and it considered just a limited number of ionic 
rotational states. In this study, we use accurate vibrational wave functions and a larger number of 
possible rotational states of the ground vibrational level. The DR rate coefficient obtained is 
found to agree better with the experimental data from storage-ring experiments than the previous 
theoretical calculation. We present evidence that excited rotational states could be playing an 
important role in those experiments for collision energies above 10 meV. The DR rate coefficients 
calculated separately for ortho- and para-Hjj~ are predicted to differ significantly at low energy, a 
result consistent with a recent experiment. We also present DR rate coefficients for vibrationally- 
excited initial states of Fhj~, which are found to be somewhat larger than the rate coefficient for the 
ground vibrational level. 

PACS numbers: 34.80.Ht 34.80.Kw 34.80.Lx 
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I. INTRODUCTION 



Dissociative recombination (DR) of the simplest polyatomic ion Hjj 



H+ + e 



H 2 + HorH + H + H 



has been studied for several decades both in experiment and theory 
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(i) 

The measured 



12J for electron energies below 1 



rate of the reaction is relatively fast 
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eV, which was eventually attributed to the strong Jahn- Teller coupling between vibrational 
motion of the ion and the incident p-wave in thej^lectronic continuum 3), 0]. Current 
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rere is general agreement between theory 
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13] and most recent experiments 



as 



12 for the DR rate coefficient in HT. However, the detailed energy dependence 
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of the theoretical rate coefficient |4|, [5| in the range 0-2 eV exhibits differences from the 
rate coefficient measured in recent high resolution storage ring experiments . One 

prominent point of disagreement is the much more pronounced resonance structure in the 
theoretical rate coefficient; this plethora of resonances is associated with Rydberg states of 
the neutral molecule Hg. Although similar resonance structure is visible in the experimental 
data, it is much less pronounced. 

Dissociative recombination of is a four-body problem: Because the DR process starts 
with an electron-H^ collision, one must account for the motion of the electron and its cou- 
pling with the molecular degrees of freedom. In order to represent Jahn- Teller coupling 
between electronic and vibrational motion, one must take into account at least two degrees 
of freedom of vibrational motion (two hyperangles in our approach). Finally, the third vibra- 
tional coordinate (the hyper-radius) should be included in order to describe the dissociative 
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channel. Therefore, the theoretical approach includes several ingredients [4], [5| and a number 
of approximations were made to simplify the calculations in our previous theoretical study 

One of the important approximations used in Refs. % |5) is the adiabatic hyperspherical 
approximation: motion in the hyper-radial coordinate was treated as adiabatic compared to 
motion in the hyperangles, i.e. motion in the hyperangles was considered to be much faster 
than motion in the hyper-radius. Correspondingly, vibrational states Q v of the ion have 
been represented as simple products of hyperangular and hyper-radial functions. Somewhat 
surprisingly, this approximation describes reasonably well (at about the 1% level) vibrational 
energies of the ion for the low vibrational states in low-lying hyperspherical potential curves 



(see Table H]). For excited vibrational states with energies higher than 1 eV above the 
ground rovibrational level, the error increased to the vicinity of 10 meV. Such highly excited 
vibrational states may support Rydberg states of the neutral H3 molecule in the energy region 
corresponding to the low-energy DR process under consideration. Accordingly, the absolute 
positions of such Rydberg states were presumably calculated with an error of about 10 meV. 
We should also mention that the highest vibrational levels that must be included in our 
theoretical treatment (in order to represent the initial dissociation of the neutral molecule) 
have around 4 eV of vibrational energy. Another source of error in the position of Rydberg 
states is the neglected energy-dependence of quantum defects used in Refs. 0,0]. Analysis 



of an ah initio calculation 



14j shows that the quantum defects are in fact only weakly energy- 



dependent close to the equilibrium geometry of the molecular ion. The maximum error in 
positions of Rydberg states for n « 2 — 3 associated with this approximation is estimated 
to be around 10-15 meV, at least for Rydberg states close to that equilibrium configuration. 
Certain Rydberg states appear in the DR spectrum as resonances and, therefore, play an 
important role in the detailed comparison of the theoretical DR rate coefficient with the 
high-resolution experimental data. On the other hand, the DR rate coefficient that has 
been thermally averaged over a Maxwell-Boltzmann distribution should be comparatively 
insensitive to the detailed positions of Rydberg states, as long as the average level density and 
the resonance widths are described correctly. Indeed, the theoretical thermal rate coefficient 
calculated in the previous study 4j agrees with the experiment ally- measured thermal rate 
coefficient. 

The adiabatic hyperspherical approximation implemented previously for the ionic vibra- 
tional eigenstate calculation neglected all non-adiabatic coupling between different adiabatic 
channels. This approximation could also adversely affect the calculated DR rate coefficient 
since the main DR mechanism in Hjj~ is indirect: neglect of some non-adiabatic effects might 
therefore be expected to cause an underestimation of the theoretical DR rate coefficient. 
However, the most important coupling responsible for the high DR rate coefficient of Hjj~ is 
non-Born-Oppenheimer Jahn- Teller coupling, which has been accounted for in the previous 
study. 

Finally, in the previous study only a few rotational levels of the initial state of the ion 
were included up to N + = 3. (The rotational state (3,1) was not included. Here and below 

n 

we use the generally accepted notations for rotational states of the ion [15J, where the two 
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numbers in () are the total angular momentum N + of the ion and its projection K + on the 
ionic symmetry axis.) The inclusion of a larger number of rotational states has since been 
recognized as possibly important due to the following: If the electron energy is high enough, 
there is a large probability that the ion is excited into a higher rotational level {4, 16]. 
In fact, the probability of this rotational excitation is high enough to be competitive with 
the DR process. Thus, for electron energies above 10 meV, higher rotational states of the 
ion could be populated in the storage ring even though, initially, the ionic beam has been 
prepared in the ground rovibrational state. Excitations can hap pen, for instance, when the 
ions pass through the toroidal region or the electron cooler 10 ] . 

The present study improves upon two of the three approximations discussed above. First, 
we calculate accurate vibrational states of the ion instead of relying on the adiabatic hy- 
perspherical approximation; the adiabatic hyperspherical representation is still utilized, but 
since channel coupling is now included, the calculated eigenspectrum can be made arbitrarily 
accurate, in principle. Second, we account for a larger number of initial rotational states. Fi- 
nally, we calculate the rate coefficients of DR processes that start from two different excited 
vibrational states of FFj\ 

The article is organized as follows. Section II briefly summarizes our theoretical approach 
and, in particular, the method to obtain accurate vibrational energies and wave functions 
of the ion. Section II I II discusses the various averaging procedures that must be performed 
on the raw theoretical DR rate coefficient in order to compare theory with the data from 
existing storage ring experiments. Section [TV] presents our results and Section [V] gives our 
conclusions. 



II. THEORETICAL APPROACH 

The procedure of calculation of the cross-section is very similar to that described in Ref. 
[n]]. Here we briefly summarize the main steps. 

The treatment is based on construction of the multi-channel electron-H^ scattering matrix 
S. The asymptotic channels \i) of the matrix corresponds to different rovibrational levels of 
the ion. After a collision with the electron the rovibrational state of the ion may change: 

e - + H 3 + W^ e -+H 3 + (0. (2) 
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However, the collision conserves the overall symmetry of the system. (That is, S is diagonal 
in the irreducible representation T to t of the symmetry group D^h of the Hamiltonian.) We 
assume that the electronic state of the ion is the ground 1 A[ state. Since protons are 
fermions, for Hjj" in the ground electronic state, the allowed irreducible representations of 
the total nuclear wave function (including space and nuclear spin coordinates) are A 2 and A 2 
of the D^h symmetry group. When the scattering matrix is constructed, highly excited 
vibrational levels of are included. These vibrational levels in fact represent the discretized 
vibrational continuum, and they have finite lifetimes with respect to dissociation. The energy 
of the incident electron is not high enough for these levels to represent open channels for 
dissociative ionization, but the electron can be captured into the Rydberg states attached to 
the highly excited levels. If this happens, the system dissociates instead of ionizing, because 
such low principle quantum number Rydberg states are locally open for dissociation. This 
causes the electron-ion scattering matrix S to be non-unitary, and the 'defect' from unitarity 
of the relevant columns of S can be identified with the dissociation probability. 

The total wave function of the ion-electron system is constructed by taking into account 



all symmetry restrictions determined by the 
system, which is discussed in detail in Refs. 



;wo allowed irreducible representations of the 
4|, |5|. The wave functions of the A' 2 and A' 2 ' 



irreducible representations of the channel states \i) are constructed from the product of the 
rotational, vibrational, nuclear, and electronic degrees of freedom of the system (see also 
Eqs. (2) of Ref. [51): 

\l) = P^total , 
§total = $rot®vib®ns$el • (3) 

The operator P s y m projects the product ^totai on the corresponding irreducible represen- 
tation, A 2 or A 2 . Each factor in the second equation is calculated by the diagonalization 
of the respective Hamiltonian, except $ e ;. In the above equation, $ rot is the rotational 
wave function of the ion, Q v n, is the vibrational wave function, Q ns is the nuclear spin wave 
function, and $ e ; represents the wave function of the incoming electron. The functions 
$ e z do not diagonalize completely the clamped-nucleus electronic Hamiltonian; the matrix 



of the electronic Hamiltonian represented 
Jahn- Teller coupling in H3 [ij], [b 
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jy $ e ; has non-diagonal elements responsible for 



2l|. 



The rotational, nuclear spin, and electronic functions in the product Qtotai are represented 



in the same way as in our previous study, Ref. 4J. But the vibrational part is obtained 



differently here, because we no longer utilize the hyperspherical adiabatic approximation 
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24i | . In that approximation, the non-adiabatic couplings between different hyperspherical 



adiabatic states of the same vibrational symmetry were entirely neglected, an approximation 
that works reasonably well for low vibrational levels of Hjj" ion. In this study we use much 
better vibrational states calculated using the slow variable discretization (SVD) approach 
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261 ]. where non-adiabatic couplings are taken into account. Table [J compares accuracy 



of calculations of the present approach with the hypers phe rical adiabatic approach. The 



calculations of H^" vibrational states by Jaquet et al. 



Oppenheimer ionic surface that we use here 
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271 ]. obtained for the same Born- 



are taken as an 'exact' reference. 



Once ^totai are known, the matrix elements 5» »/ are calculated from the scattering matrix 
5a,A'(Q) depending on three distances between protons in H3 and describing the e~+H^ 
collision in the molecular frame, where the appropriate quantum numbers are projection A 
of the electronic orbital momentum on the ionic principal axis and the set of internuclear 
coordinates (Q). As in the previous study {4], 5], we consider only the l p— wave' of the 
electron, when it moves beyond the range of the ionic core. The nonspherical nature of the 
electron-ion interaction potential undoubtedly mixes other electronic orbital momenta when 
the electron is at short range, but scattering calculations have shown that the probability 
for an incident p — wave electron to scatter into an s or d orbital momenta is quite low 
in this near-threshold energy range. The matrix Sa.a'(2) i s obtained from the reaction 
matrix K\ \/{Q) given by formulas in Refs. 18L Il9j]. As was mentioned above, the electronic 
Hamiltonian and, correspondingly, the matrices related to body-frame scattering, i.e. S 
and K, are not diagonal. The nonzero off-diagonal elements K^__ x /Sx-\ are due to the 
Jahn- Teller coupling. 

When the total energy of the e~+ELt system is not high enough for all the channels |z) to 
be energetically open for ionization, the usual situation, the physically meaningful scattering 
matrix S phys (E) is obtained from by the standard closed-channel elimination procedure 
of multi-channel quantum defect theory (MQDT) (29), 30]. The dissociative recombination 
rate coefficient is then calculated using the unitarity 'defect' of the corresponding columns of 
S phys 0, 5I . In order to compare our results with the data from storage ring experiments, we 
carry out a number of averaging procedures in order to model the experimental conditions, 
as detailed below. 



6 



h ■ 
V\V2 i irrep. 


adiab. approx. 


Dvu caic. 


Jaquet et al 


27 


1 


nnO a 

UU Ji\ 


U 


n 
U 





1U J±\ 


OlOO 


oil / .0 


3178.15 


Uz J±i 


A 7^ A 
4 1 04 


A 777 Q 

4/ m .y 


4778.01 


ZU /ii 


UZ 1 


uzuu.u 


6261.81 


noO /I 
Uo /ii 


7"?S9 
f OOZ 


707c 
/ Z 1 0.0 


7285.32 


i oO /I 
Iz /li 


f 04o 


7779 Q 

( ( ( z.y 


7769.06 


U4 ^il 


oa It) 


oyyo.o 


9000.58 


QnO a 
oU 


»z4o 


yzoo.o 


9252.08 


lO 74.1 


1 ni 9Q 
luizy 


yyoo.u 


9963.98 


99° /I, 

ZZ s\\ 


1 n/L9n 

lU^iZU 


iuoyo. 


10590.51 


UO l\\ 




1 noi 9 
iuy iz. 


10915.47 


nil?? 

Ul X!/ 


9^1 R 
zOlO 


9^91 1 

ZOZl.l 


2521.20 


H92 I? 
UZ i!/ 


0UU1 


Anna a 

4yyo.o 


4997.73 


1 1 1 TP 


0004 


oooz.y 


5553.95 


n^l TP 
UO -C/ 


uy 1 


KQQQ 9 

uyyy.z 


7005.81 


1 02 771 

lz £/ 


78Q7 

1 oy / 


1 oOO.Z 


7869.82 


91 1 TP 
zl 


QA 7C 


o4o 1 .0 


8487.53 


n/|2 I? 
U4 III 


yioi 


yuyu.o 


9112.90 


lO Hi 


y 1 00 


yu4y .z 


9653.42 


U4 x!/ 


youz 


noon 9 
yyyy.z 


9996.72 


992 TP 
ZZ 


1 nfi77 

1UU ( f 


1UU4U. 


10644.59 


nc^l TP 

Uo 


1091b 


1082Y. 


10862.46 


31 1 TP 

Ol J_y 


1 1 9fi c i 
1 1ZUO 


1 i04y. 


11322.31 


1 4 2 jp 
l^i 7j 


1 1 73Q 
11 / oy 


1 1UOD. 


11657.69 


05 5 £ 




12078. 


12078.43 


03 3 A 2 


7482 


7493.2 


7491.89 


13 3 A 2 


10243 


10209.7 


10209.55 



TABLE I: Accuracy test of the adiabatic hyper spherical approximation and the improved coupled- 
channels hyperspherical calculation adopted for the computations presented in this paper. Specif- 
ically, this table compares several vibrational energies in cm" 1 calculated in the present approach 
with the older adiabatic approximation results and those taken from a full three-dimensional diag- 



III. CALCULATION OF THE RAW DR RATE COEFFICIENT AND ITS AVER- 
AGE 



The rate coefficient for dissociative recombination depends on the initial rovibrational 
state \rv) of the ion. Using the defect of unitarity of the physical scattering matrix S phys , 
the DR rate coefficient a rv for a particular rovibrational state \rv) is given by {4] 



where E e i is the kinetic energy of the electron at infinity; the channel index i' at S corre- 
sponds to the initial \rv) state. The scattering matrix S is calculated separately for each 
total angular momentum N of the ion-electron system, N + is the angular momentum of the 
initial \rv) state. Notice that several values of N may contribute to the DR rate coefficient 
for the given state \rv). 

We now address the way we account for the experimental conditions in storage ring 
experiments, especially the experimental distribution over relative velocities of the ion and 
electron. In the storage ring experiments, the distribution is not uniform: the parallel 
component u\\ of the e~— Hjj~ relative velocity u has a smaller distribution width than the 
perpendicular component u±. The rate coefficient a rv (E e i) also depends on the rotational 
level. In the experiments, the initial vibrational state is usually the ground state {00°}, 
but several rotational levels are typically populated. The population of different rotational 
levels is accounted for by introducing a finite rotational temperature T rv of hht , though it 
should be remembered that this assumption that the ions are in thermodynamic equilibrium 
at some T has not been explicitly confirmed experimentally. Present generation storage ring 
experiments measure the DR rate coefficient as a function of the parallel component En 
of the total relative energy E e i of the ion and electron. The average over the non-uniform 
electron velocity distribution is then given by the following formula [17J 
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El^f 1 - £ ^ s {E el )^{E el )\ 

N \ i=l,N a J 



(4) 



a sr (E\\) = — ^ / du\\ dE±a rv [(v\\ + u\\) 2 /2 + E_l\ w sr (rv,T rv ), (5) 




where the normalization constant N sr and the statistical factor w sr (rv,T rv ) are 



w sr (rv,T rv ) = (21 + 1){2N + + 1) exp 
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In the above equations, AE± and AEu are distribution widths (measured in energy units) 
for the parallel u\\ and perpendicular u± components of the relative velocity v — m + it\\ + 
u±; v\\ represents the center of the velocity distribution, i.e. velocity at which the actual 
measurements are made in the storage ring experiments: E\\ = ujf/2. The perpendicular 
component of the energy is E± = u\/2. I is the total nuclear spin, which can be | or 
| depending on the rovibrational state \rv); E rv is the energy of the Hjj" rotational state 
(assuming that the vibrational state is always the same). The sums in the above equations 
are over all possible rovibrational states \rv) including all symmetries of the rovibrational 
states that can be populated at a given rotational temperature T rv . 

To compare with the storage ring experiments, one also must take into account the so- 
called toroidal effect, which is due to the geometry of the merged electron and ion beams, as 
there are two regions where the electrons are bent into or out of a trajectory that is parallel 
to the ions. The experimentally-observed rate coefficient with the toroidal effect correction 
is [17|: 

2 phend 

a t or{E\\) = a sr (E\\) + — a sr {E\\(x))dx . (7) 



The function E\\ (x) and the length hend account for the geometry of the merged electron and 



ion beams. A detailed discussion can be found in Ref. 



17| . After the averaging procedures 



described above, the DR rate coefficient atp r (E\\) can be compared with the raw experimental 



data from the storage ring experiments |9j, [10|, [12 ] . 



The thermally averaged DR rate coefficient relevant to a situation in which the ions 
and electrons are in common thermal equilibrium at temperature T is calculated from the 
following integral: 

a th (kT) = — a rv (E e j)w(rv, kT) y/E^dEg (8) 

Nth Jo ^ 

where the normalization constant N t h and the statistical factor wth{rv, kT) are 

POO 

N th = / ViBttKiT)^^, 
Jo rv 

w th {rv, T) = {21 + 1)(2N + + i) e -E rv (R)/kT e -E el /kT _ (g) 
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FIG. 1: Comparison of experimental j], [icl [l^ | (black circles and red diamonds) and present 
theoretical (solid line) dissociative recombination rate coefficients. In the theoretical calculation, 
the rotational temperature is T rv =1000K, the widths A.-E'ii of the parallel component of the electron 
energy is 0.1 meV, and AE± =2 meV. 



IV. RESULTS 



Figures [U [2J, [3], and 0] summarize the computational results of the present study. Fig. 
compares the experimental DR rate coefficient from a recent storage ring experiment 
J, [l(3] with the present calculation. Overall agreement with experiment is better than in 
the previous theory (see Fig. 6 in Ref. 13]). This is due to two factors: In the present 
treatment, we use more accurate vibrational wave functions, which are calculated using 
SVD, i.e. without the hyperspherical adiabatic approximation. The second improvement 
is due to the larger rotational temperature T rv and the larger number of rotational states 
that are taken into account in the averaging formula Eq. (jSJ). The energy of the highest 
rotational level (5,1) included in the present calculation is 1250.3 cm -1 15j. The energy 
of the lowest state (1,1) allowed for Hjj~ is 64.1 cm -1 above the symmetry-forbidden state 
(0,0). In Table [Til we show the partial contributions of the rotational states with different 
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N + to the total DR rate coefficient with T r „=300 K for four different energies 0.001, 0.01, 
0.1 and 1 eV. Although the relative population of the N + = 5 rotational states is about 3% 
at T rv =300 K, there can be accidental cases, at certain energies where other contributions 
happen to be small, where it can become an important contributor to the observed DR rate. 
For example, at a collision energy 0.0997 eV, the cumulative contribution of the states with 
N + = 5 contributes 14% of the calculated DR rate. 
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FIG. 2: This figure presents calculations with different rotational temperatures in Eq.©. The 
theoretical DR rate coefficients obtained at higher rotational temperatures, e.g. 300K-1000K, 
agree better with both the CRYRING and TSR experiments. 



Figure [2] demonstrates the theoretical DR rate coefficients we obtain for different ionic 
rotational temperatures T rv . The results with higher T rv agree better with the experiments 
than when the experimentally estimated rotational temperature is adopted. While this 
better agreement at a higher rotational temperature could be fortuitous, our results suggest 
that it is worth exploring whether the rotational temperature in both of the recent storage- 



ring experiments [9, 



10 



12| might be larger than 40 K or 13 K, respectively. (40 K and 13 K 



are the estimated rotational temperatures in the two experiments). This conclusion would 



conflict with another suggestion by Kreckel et al. 



12J, that only the two lowest rotational 
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Energy (eV) 


Total DR rate coefficient 

(cm 3 /s) 


N+ = l 


N+ = 2 


N + = 3 


AT+ = 4 


N + = 5 


0.00101 


2.19 x 10~ 7 


0.295 


0.198 


0.364 


0.126 


0.017 


0.0103 


4.79 x 10~ 8 


0.388 


0.139 


0.301 


0.136 


0.035 


0.0997 


6.30 x 10~ 9 


0.157 


0.238 


0.266 


0.196 


0.143 


1.02 


7.10 x 10~ n 


0.324 


0.201 


0.315 


0.111 


0.049 



TABLE II: Partial fractional contributions to the DR rate coefficient from individual ionic an- 
gular momenta, N + = 1, • • • ,5, at four different energies. The DR rate coefficient is calculated 
for T r „=300K. In the table, the DR rate coefficient and partial contributions are calculated with- 
out the toroidal averaging. For larger temperatures, higher rotational states can have important 
contributions at certain electron energies. 

states are populated in the recent storage ring experiments fly , [12[ and that, therefore, 
the rotational temperature T rv is about 13 K — 40 K. In fact, it was demonstrated previously, 
that if the electron energy is high enough, the electron-ion collision might not only cause DR 
or result in an elastic collision, but it can also result in rotational excitation of the ions when 
they circulate in the storage ring 0, [if]]. In all the tests we have carried out for temperatures 
in the range 13 K-40 K, the calculated DR rate coefficient has pronounced structure due 
to Rydberg states present in the raw rate coefficient of Eq. HI The experimental DR rate 
coefficient has some structure, but it is less pronounced than our averaged and convolved 
theoretical rate coefficient. Since the resonances due to the Rydberg states are smeared 
out in the experiments, it suggests the possibility that in the experiment there could be 
an additional source of broadening. The broadening could arise from a higher rotational 
temperature, from a broadened electron energy distribution, from additional broadening 
associated with the toroidal region, or perhaps from something else. 

An alternative possibility that cannot be ruled out is that our theoretical treatment 
might have underestimated the resonance widths. For the resonances that dominate the 
DR rate, the predissociation partial width is larger than the autoionization partial width, 
and under those conditions, the calculated DR rate is comparatively insensitive to changes 
in the predissociation linewidth. Thus, it would be a valuable benchmark for experiments 
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(or other, improved theories) to determine the predissociation partial widths of individual 
resonances above the ionization threshold, to provide a direct test of the accuracy of our 
present calculations at the level of spectroscopic accuracy. 

One possible source of rotational excitation could be the repeated circulation of the 
molecular ions through electron cooler during the ramping of the cathode voltage [l(3|. The 
authors of Ref. [lO] deduced that the rotational temperature is 40 K based on theoretical 
cross-sections of the rotational excitation of Hj~ given in Ref. 3JJ. Since then, the cross- 



sections have been reconsidered and corrected 16j. Correspondingly, we revisit here the 
arguments of Ref. 1(| based on the new inelastic probabilities determined by Ref. fl6 l. 
If we take the (1,1)— >(2,1) rotational excitation cross-section to be 710 A 2 from Ref. 16j 
instead of 210 A 2 from Ref. [31], we obtain the relative population 7.5% instead of 2.2% in 



Ref. 31]. If we take the (1,0)— >(3,0) rotational excitation cross-section to be 270 A 2 from 
Ref. 16] instead of 120 A 2 from Ref. [3lJ, we obtain the relative population 4.6% instead 
of 2.0% in Ref. [3j|. The relative population 7.5% of the (2,1) states corresponds to the 
temperature T rv =96 K, the relative population 4.6% of the (3,0) states gives T rv =210 K. 
Both values of the temperature are significantly higher than the values (40 K and 13 K) 
quoted in the experimental papers, although not as high as 1000K, in our present estimation. 

The disagreement between the current theory and experimental data around energy 
E\\ = 6 meV (see Fig. [Q could be caused by errors in the calculated positions of Ryd- 
berg states in that region, which are attached to highly-excited rovibrational levels of the 
ion. In our calculation, we use energy-independent quantum defects even though they de- 
pend weakly on the principal quantum number. In addition, the accuracy of calculation 
of energies for the highly-excited rovibrational levels may be of the order of 6 meV. Note 
that there is also a disagreement with experiment in the region of very small energies, below 
0.2 meV. This region is well below the net effective energy resolution of the experiment 
(AE± =2 meV). The experimental rate coefficient appears to behave there as E\~ l l 2 , which 
is the same total energy dependence expected for the raw, unconvolved DR rate coefficient. 
However, in our theoretical calculation, we of course included the convolution according to 
the perpendicular energy distribution with the width AE± =2 meV. This makes the con- 
volved theoretical DR rate coefficient become essentially flat at very low energy E\\ <C 2 
meV even though the raw theoretical rate coefficient also grows as E- 1 ' 2 at very low energy. 
This suggests that in the experiment, the distribution of relative electron energies could 
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be even more complicated than discussed above, and in particular, the resolution might 
be even better at very low energies than the quoted energy resolution. Another possible 
explanation for the disagreement is that a Rydberg resonance exists in H 3 at an energy just 
above (+0.3 meV) the (1,1) rotational state of the ion. The predissociation linewidth of the 
resonance must also be of the order 0.3 meV. We have made a simple test calculation in 
which we artificially tuned one of the Rydberg resonance to be placed just above the (1,1) 
ionization threshold. The resulting theoretical DR rate coefficient looks very similar to the 
experimental DR rate coefficients at energies below 0.5 meV. Therefore, the sharp increase 
of the experimental DR rate coefficient for energies below 0.3 meV requires additional con- 
sideration. Similar discrepancies have been observed between DR theory and experiment at 
very low parallel energies, in other systems such as LiH" 1- ,!^] so this could be a systematic 
issue for theory and experiment to confront which extends beyond the Hjj" system alone. 

Figure [3] compares the present theoretical DR rate coefficient (dashed grey curve) with 
the recent TSR storage ring experiment by Kreckel et al. 12j. In the TSR experiment, the 
parallel and perpendicular energy resolution parameters A- Ey =25/xeV and AE±=0.5 meV 
are slightly smaller than in the CRYRING experiment 0, [lo|. Therefore, the theoretical 
curve shown in the figure has been correspondingly convolved using these parameters. These 
calculations have also assumed that the target ion rotational temperature T rv = is equal to 
1000 K. The overall agreement between theory and experiment is good except in the energy 
region below 0.15 meV already discussed above. As one can see, when the width AE± is 
decreased from 2 meV in Fig. [I] to 0.5 meV in Fig. [3j the theoretical rate coefficient agrees 
better with the sharp increase of the experimental DR rate coefficient for energies below 0.5 
meV. 

In the TSR experiment, Kreckel et al. 12] . observed the dependence of the DR rate 



coefficient on the nuclear spin of Hg~. They found that para-Hjj" has a larger DR rate 
coefficient than ortho-H^ for low energies (< 0.5 meV). The previous theory {4] has predicted 
different DR rate coefficient for para-H^" and ortho-hLt: At low energies the theoretical 
DR rate coefficient for ortho-H^ was larger than for para-H^", i.e. opposite to what was 
observed in the experiment by Kreckel et al. We now revisit this issue in the context of our 
new and presumably improved theoretical description. Figure [3] shows the separate ortho- 
H3" and para-H^ DR rate coefficients calculated in the present treatment. The para-H^ 
rate coefficient is significantly higher than the rate coefficient obtained for ortho-ILt. This 
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FIG. 3: This figure compares the theoretical DR rate coefficient to the high-resolution storage 
ring experiment of Kreckel et al. [121 ] carried out at TSR. The experimental resolution parameters 
are A-Eii and AE± are 25/ieV and 0.5 meV respectively. The theoretical curve shown has been 
calculated with these parameters and rotational temperature T rv =1000 K. The figure also shows 
the theoretical DR rate coefficients calculated separately for ortho- and para- configurations of 
with the same parameters A-Em, AE±, and T rv . 



dramatic difference between the ortho-H^ and para-H^" rate coefficients obtained at low 
electron energies and those of our previous theoretical study appears to result from slightly 
different positions of the calculated Rydberg H3 states whose energies lie close to the (1,1) 
and (1,0) ionic rotational states. 

Figure H] presents our theoretical DR rate coefficients obtained for a target ion that 
is initially in an excited vibrational state. These calculations have been carried out for 
the first {01 1 } and second {10°} excited ionic vibrational states. (In { | we specify the 
vibrational quantum numbers of the ion using the normal mode notation 15j.) The energy 
of the lowest rotational state (0,0) for {01 1 } is 2521.4 cm -1 , and the energy of the lowest 
rotational state (1,1) for {10°} is 3240.7 cm -1 15j. The DR rate coefficient for the two 
excited vibrational levels of the ion is higher than the DR rate coefficient for the ground 
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state, which is reasonable, considering that a similar qualitative increase of the DR rate 
coefficient was previously observed in both theory and experiment for diatomic ions. 




Electron energy, E (eV) 

FIG. 4: Comparison between the theoretical DR rate coefficients for the H^j" ion prepared in the 
ground and excited vibrational levels. At low energy, the DR rate coefficient for the vibrationally 
excited ion is significantly larger than for the ion in the ground state. This result is consistent with 
trends observed for DR rate coefficients in diatomic ions, where the DR rate coefficient typically 
increases with vibrational excitation. In the legend, the numbers in parentheses have the same 
meaning as in Fig. [2j 

Finally, Fig. [5] compares the theoretical and experimental thermal rate coefficients. The 
theoretical rate coefficients are obtained directly from the raw theoretical data using Eqs. 
(JHJ) and (jUj). Thus, the toroidal effect and the finite widths AE±, A£j|, kT rv are not present 
in these theoretical results. In fact, our calculation shows that the inclusion of the toroidal 
correction increases the thermal rate coefficient by about 20 % approximately uniformly for 
all energies. Thus, the agreement with experiment is good. The theoretical thermal rate 
coefficient at 300K is 5.6 x 1CT 8 cm 3 /s. 
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FIG. 5: The present theoretical thermal rate coefficient for dissociative recombination of Hg is 
compared with th e experimental rate coefficient deduced from the storage ring experiment of 
McCall et al. 
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V. CONCLUSIONS 



In summary, we would like to emphasize the following results from the present study: 
We have calculated the rate coefficient of H^" dissociative recombination using an im- 
proved description of the ionic vibrational states, and including more target rotational 
states. The resulting theoretical rate coefficient agrees reasonably well with two recent 



storage ring experiments 



121 ] . The agreement with experiment has been improved 
over that achieved in the previous theoretical study |4j, |5| ; in particular, the present results 
may point to a resolution of the largest previous discrepancy, in the energy range from 
0.04 eV to 0.15 eV. However, the improved agreement with the experimental data in that 
energy range was only obtained when we assume that the rotational temperature T rv of 
it is significantly larger than 40 K and 13 K, the values given in the experimental study 
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121 ] . Since no direct measurement of the temperature has actually been made inside 



the storage ring in those experiments (except indirectly for zero-energy collisions), there is 
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a possibility that the ions get rotationally excited before the DR rate coefficient measure- 
ments are conducted. It was shown previously that the probability of rotational excitation 
of the ion by electrons is comparable to or larger than the DR probability, at energies where 
rotational excitation is energetically allowed. Thus, the rotational temperature in the ex- 
periment could be larger than 40 K or the ions might not even be in thermal equilibrium at 
any temperature. Thus, it would be desirable to monitor the rotational temperature during 
the DR measurement. Another possibile way to explore this effect is to artificially increase 
the temperature of the electron cooler or the width AE± during the DR measurements and 
ascertain the temperature at which the DR rate coefficient starts to become sensitive to the 
temperature. 

We have calculated the DR rate coefficients for separate ortho- and para-configurations of 
H^ . At energies below 10 meV, the DR rate coefficient for para-Hjj" is an order of magnitude 
larger than for ortho-H^ . The experiment also shows that the para-Hjj" DR rate coefficient 
is larger. However, since the ortho-/para-ratio in the experiment is not known it is not clear 
what is experimental DR rate coefficients for pure para-Hg and ortho-Hg . Our previous 
calculations stressed that the ortho-para ratio of DR rates at very low collision energies 
should be used with some caution, because the rates at energies below 100K begin to get 
very sensitive to the specific resonance positions at the meV level. Those cautionary remarks 
are still applicable to the present results for the ratio of ortho and para DR rates at low 
energy. However, if the present order-of-magnitude difference of the low energy DR rate 
survives future improvements in theory and experiment, it will be interesting to explore 
possible implications of this difference for the chemistry of interstellar clouds. 

Finally, the calculated DR rate coefficient for ions prepared in excited vibrational states 
is larger than in the ground vibrational state. Our new theoretical values for the DR rate 
coefficients of excited vibrational states will hopefully be tested one day in a storage-ring 
experiment. Currently, experimental DR measurements are made after the ions have been 
cooled. In principle, it seems possible to carry out this measurement using vibrationally-hot 
ions and, to monitor the DR rate coefficient as a function of the vibrational temperature. 
Such an experiment might give deeper insights into the energetics and target state de- 
pendence of the DR process, and it could then be compared with the present theoretical 
predictions. 
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